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Physical mechanisms that can influence rotation curves of spiral galax- 
ies are discussed. For dark matter studies, possible contributions due to 
magnetic fields and non-Newtonian gravitational accelerations should be 
carefully acconted for. We point out that magnetic fields are particularly 
important in outermost parts of the disk. In the framework of general rela- 
tivity the physical reason of an enhanced gravity in spiral galaxies depends 
on the assumed metric. The additional gravity is provided for Schwarzschild 
metric by nonluminous mass, whereas for Vaidya metric [1] by emission of 
radiative energy. In the latter case the non-Newtonian acceleration displays 
l/r behaviour. Also matter flows contribute to non-Newtonian gravity. 

PACS numbers: 98.52.Nr.,98.62.Dm. 



1. Introduction 

The best evidence of enhanced gravity in galaxies is provided by flat 
rotation curves of spiral galaxies, which do not decay in a Keplerian way even 
far from the rotation axis. From simple Newtonian formula for centripetal 
acceleration, 



■4 ot _ GM{r) 



(1.1) 



one finds that the total mass within radius r grows with r as M(r) ~ 
rVr t- The linear growth of mass is customerly attributed to an invisible 
component, referred to as dark matter. 
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The problem of dark matter has its beginning in the observational de- 
termination by Zwicky of dynamical mass of the Coma cluster of galaxies. 
The gravitational mass inferred by Zwicky from the motion of individual 
galaxies in the cluster exceeded by a factor of a few hundred the mass ob- 
tained by measuring luminosities assuming typical value of mass to light 
ratio. Later a discrepancy between dynamical and luminous mass has been 
found in spiral galaxies and galaxy clusters. 

Our aim here is to point out that the conclusion as to the existence of 
dark matter inferred from rotation curves is not inescapable, but based on 
some unspelled assumptions. It holds in Newtonian gravity provided any 
role of magnetic fields is negligible. In the framework of Einstein's gravity 
for it to hold one implicite assumes space-time geometry to be given by the 
Schwarzschild metric. It is often assumed that galactic gravitational field, 
as very weak one, is adequately described by Newtonian gravity. We will 
show, employing Vaidya metric, that the inverse problem, of reconstruct- 
ing galactic gravity given rotation velocity, has also other solutions. One 
encounters here ambiguity, which can only be resolved by physical input. 

Recent observations of dearth of dark matter in elliptical galaxies [2] sug- 
gest that there may be more unknowns involved in this problem. One should 
also consider non-gravitational origin of the above discrepancy, namely due 
to magnetic fields in galaxies. The role of magnetic fields is likely very im- 
portant in the outer disk region, where the galaxy rotation is detected by 
tracing hydrogen clouds. 

In order to firmly infer the amount of dark matter in spiral galaxies one 
should subtract contributions to rotation curves generated by other forces 
(i.e. magnetic fields) and processes. It is certain that such contributions 
exist as in many galaxies rotation curves show wiggly structure, as e.g. in 
our Galaxy. Such a structure cannot be produced by WIMP gravity, as 
density of WIMPs is a monotonically decreasing function of distance from 
the center of the distribution. Magnetic influence and energy-flow-generated 
gravity can easily account for undulations of rotations curves. However, 
immediately arises a question how much such forces/processes contribute 
to the bulk of rotation curves. Before gauging this influence the inferred 
amount of dark matter in a galaxy is subject to substantial uncertainty. 

In the next section the problem of flat rotation curves of spiral galaxies is 
briefly reminded. In sect. 3 the role of magnetic fields is discussed. In sect. 4 
gravity due to radiation flow is discussed with the use of Vaidya metric [1]. 
Finally in the last section we summarize important points once again. 
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2. Enhanced gravity in galaxies 

One can formulate the problem of flat rotational curves precisely as fol- 
lows: there is too much gravity compared to mass we can account for by 
counting stars and measuring the amount of gas in galaxies. It is a nonrela- 
tivistic custom to attribute this enhanced gravity to invisible and (almost) 
undetectable matter. In a popular Cold Dark Matter model, invisible mass 
is due to hypothetical weakly interacting massive particles - WIMPs. 

In general relativity, which is supposed to be the theory of gravity, not 
only mass is capable of generating gravitational field, but also energy or 
radiation flows induce gravity. Interpreting gravitational acceleration in 
Newtonian terms, a = —GM(r)/r 2 r, for spherical symmetry one implicite 
assumes the Schwarzschild interior metric 

ds 1 = e u dt 2 - e x dr 2 - rW, (2.1) 

with e~ x = e v = 1 — 2M(r)/r for matter with negligible pressure, described 
by a dust equation of state. In the weak field limit, which is appropriate for 
galactic fields, one obtains then the Newtonian acceleration. 

Typical rotational velocities of spiral galaxies, in the flat regime, are 
of order of 100 km/s. One can thus infer the mass within radius r to be 
M(r) = 2.32 x 10 9 (v c /W0km/s) 2 r /kpcM Q . For Milky Way galaxy this 
gives within SOkpc the mass 

M M w = 3.37 x 10 n M Q for v c = 220km/s. 
This high value of mass is thought to show us that the main component of 
mass in our Galaxy, and in other galaxies, is nonluminous. Astronomers 
tried hard to detect known nonluminous astrophysical objects that could 
form an invisible population providing the missing mass. All attempts to 
account for it by dim stars, dead stars, plasma or other forms of baryon 
matter have failed. The only viable candidate at the moment is particle 
dark matter composed of WIMPs (or axions), forming an extended halo 
around galaxies. The radius of this halo is presently unknown, but some 
observations suggest it is of order of 200 kpc. Also, some cold gas in the 
form of molecular hydrogen, can contribute to dark matter, as it is very 
difficult to detect this component. 

One can briefly summarize that the dark matter hypothesis is a New- 
tonian solution of enhanced gravity problem in spiral galaxies with any 
influence of magnetic field neglected. 

3. Magnetic fields and the rotation curves 

Magnetic fields are the most common phenomenon in spiral galaxies 
where we observe fields of regular and chaotic structure. The regular struc- 
ture has azimuthal and poloidal components. The poloidal component of 
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the magnetic field is produced by the galactic dynamo effect. The azimuthal 
component is induced from poloidal component by differential rotation. 
Regular fields created by such mechanisms may reach intensities of a few to 
several hundred microgauss. 

A question arises [3,4] whether these fields have any influence on the 
galactic rotation curves. In order to answer this question we must investigate 
the Navier-Stokes equation with magnetic field, 

p[ ift + = ~ Vp ~ ^ + nM+ (e + 1 } v(v ° ^ + r ((v x §) x §) - 

(3-1) 

We assume here the gravitational field of point mass located at the cen- 
ter of galaxy with M ~ 2 • 10 44 <7, which is a good approximation for our 
exploratory calculation. In stationary galactic disk we can neglect radial ve- 
locities and viscosity and we will compare gravitational and magnetic field 
forces. Rough estimate shows that for gas clouds of density p ~ 10 -25 <?/cm 3 
at the radius of a few tens kpc, r ~ 3 • 10 22 cm, from the galactic center, and 
for magnetic fields of a few pG, B ~ 10~ 5 G, magnetic forces are comparable 
with gravitational forces. Gravitational accleration and accleration due to 
magnetic effects are, respectively, 

M ■ G o 9 , N 

— 10" 8 cms~ 2 , (3.2) 

3 • 10~ 8 cms" 2 . (3.3) 



p ■ r 

For magnetic effects to occur the gas must be partially ionized. We know 
that at least a few percent of the hydrogen in galaxies is ionized. Therefore 
we may expect, that the magnetic fields' influence on rotational curves is 
not negligible. 

The above order-of-magnitude estimate shows that the magnetic influ- 
ence is particularly important in the outermost regions of the galactic disk, 
where the density of hydrogen is the lowest. From eq.(3.3) we find that 
when density decreases by a factor of 100, magnetic fields on the scale of 
lpG can overwhelm gravity! Let us remind that most of the dark matter 
contribution to galaxy rotation comes from the outskirts. Any unaccounted 
for magnetic field contribution can completely corrupt dark mass measure- 
ment. One should also keep in mind that magnetic fields on pG order are 
expected to be ubiquitous in the intergalactic space in clusters of galaxies. 

Taking only azimuthal component of magnetic field into account and 
assuming, that it depends only on its radial coordinate we can get a simple 
analytical form of this component which will flatten the rotation curve: 
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where = v c = const is the rotational velocity of a galaxy. The solution 
of this equation reads: 



B* = +/ * 1 , (3.5) 

r 

where C is the integration constant. 

The real magnetic fields in spiral galaxies have all the components, 
poloidal and azimuthal, nonvanishing. Therefore the radial part of Navier- 
Stokes equation has the following form: 

(^) 2 _ MG 1 ( dB r dB 1 d ldB r \ 

To possible influence of such magnetic fields on rotation curves one 

should take into account the whole 3D structure of magnetic field. 
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Fig. 1. The dashed-dotted line is a solution of equation (3.8) with M = 2 • 10 44 g, 
«0 = 220- km/s, p = 10~ 25 g/cm 3 , dotted line is a solution corresponding to density 
p = 0.333 • 10~ 25 g I cm 3 . Solid curve is a function ~ l/r 



G 



VAIDYA PRINTED ON FEBRUARY 2, 2008 



4. Gravity generated by radiation flow: the Vaidya metric 

The Schwarzschild metric describes strictly speaking gravity of cold 
spherically symmetric astrophysical body, such as a planet, dead star or 
a black hole. Gravity of radiating objects such as normal stars is only ap- 
proximately described by the Schwarzschild metric. There exists, however, 
exact solution of Einstein's field equations corresponding essentially to a 
real star which emits radial flux of radiation, found by Vaidya [1]. In its 
original form, the Vaidya metric is 

ds * = > )2(1 _ ^)-1^2 _ (1 _ ^)"l dr 2 _ r 2 dn 2 (41) 

m! r r 

It corresponds to the space-time region outside the star, r > ro, where vq 
is the stellar radius, and m = m(r,t). This metric can be cast is a very 
elegant form employing the retarded time variable u = t — r, 

ds 2 = -(l-^)du 2 -2dudr + r 2 dCl 2 , (4.2) 

as shown by Vaidya in Ref.[l]. 

The energy tensor corresponding to the Vaidya metric has non-zero Tg 
component, which describes the energy outflow carried away by massless 
fields. Let us consider the energy tensor for directed flow of radiation [5] in 
the form 

T; = pv^, (4.3) 

where p is the energy density of radiation, and the fourvector is null, 
v^Vft = 0. For the radial outflow, v 2 = v 3 = 0, and T 2 2 = T 3 3 = 0. The metric 
(4.1) is a particular example of a general non-static spherically symmetric 
metric [5], 

ds 2 = e u{r ^dt 2 - e Hr ^dr 2 - r 2 dn 2 . (4.4) 
The Einstein's field equations are [5]: 



- 87 rT ° = -l + e-* ( l-^), (4.5) 



-SirTl = -l + e ^(l + -), (4.6) 



-8vrT 2 2 = -je- u (2\+\(\-v)) + ^e- x (2v"+v' 2 -\'v f +2^-^) = -8ttT 3 3 , 

(4.7) 

-87rT i = — A. (4.8) 
r 
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Let us introduce the mass function m(r, t) through e A ( r '*) = l—2m(r,t)/r. 
From the null condition, = = -e x (v 1 ) 2 +e v (v ) 2 we find e^" A )/ 2 T 1 ° + 
Tq = which gives 

e - x ' ,2 m + e~ u/2 rh = 0. (4.9) 

This allows us to express the function e v through m(r, t) and its derivatives, 
and to write the metric (4.1) in the form (4.2) given by Vaidya. 

The physical interpretation of Vaidya's metric is straighforward. In the 
weak field limit we find rh + m' = and the energy flux flowing out of a 
sphere of radius r is 



^01 



ii 

m m 



Airr 2 rh 4irr 2 
Hence m(r, t) is the radiation energy inside this sphere, 



(4.10) 



m(r,t) = J 4vrr 2 T °dr. (4.11) 



The function rh is the rate of energy emission, or total luminosity, and 
m' = 47rr 2 TQ . We should also include the radiation source, located at the 
origin, which loses energy at a rate M(0, t — r) = m(r, t). 

The most important result is the gravitational acceleration in the weak 
field limit, e v ~ 1 + v. From general expression we find for the metric (4.4) 
the Einstein's formula 

r = --u'. (4.12) 

When applied to the Schwarzschild metric in the weak field limit, with 
v = — 2M/r, it gives the Newtonian acceleration. 
For the Vaidya metric from eq.(4.6) we have 

v' = ^[(e A - 1) +e x (m'- ™)]. (4.13) 
In the weak field limit we find 

.. m' m ,aia\ 

r = — —-^2- (4.14) 

This expression shows that there appears a non-Newtonian acceleration 

m ' ^ I a ir\ 

a L = = G — , (4.15) 

r rc 

which is inversely proportional to the distance. Far from the center, 
becomes dominating. For radiating body, with energy flowing out of the 
central mass, rh < and the the acceleration (4.14) produces an attractive 
force which becomes stronger that usual Newtonian gravitation. 
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The additional gravitational attraction due to radiation emission implied 
by the Vaidya metric was first discussed by Lindquist et al.[6]. It gives 
an explicit example of non-Newtonian gravitational force resulting from 
Einstein's gravity theory for a realistic metric in the weak field limit. Thus 
the notion that Newtonian acceleration is the only weak-field limit of general 
relativity is inaccurate. 

The formula (4.14) can be generalized to the case of galactic wind which 
is a radial matter outflow, 

— • (4.16) 

Here v w i nc i is the radial velocity of the wind and m < is the mass loss 
rate due to wind. Please note that the above formula (4.15) is also valid 
for radial accretion, with radial infall velocity v r < 0. Since for accretion 
the mass increases, m > 0, the induced acceleration is also directed inward, 
as for the wind. One can thus conclude that radially oscillating shell of 
gas would always produce gravitational attraction, both in expansion and 
contraction phase. 



5. Discussion 

Astronomers tend to consider the Newtonian solution of the enhanced 
gravity problem in spiral galaxies to be the only one compatible with general 
relativity. One can encounter statements that any non-Newtonian gravita- 
tional acceleration in galaxies, as e.g. employed by Milgrom in his model 
of galaxy gravity [7], would necessarily require modifications of Einstein's 
gravity theory. We have given here an example that the statement that non- 
Newtonian gravitational acceleration, a ~ 1/r, is incompatible with general 
relativity, is not true. The acceleration (4.14) can be shown to produce flat 
rotational curves of spiral galaxies. The centripetal acceleration when the 
non-Newtonian acceleration dominates, is 

r cr 

which allows us to calculate the source luminosity L = —m. For v ro t = 
lOOkm/s L ~ W 52 erg/s. Hence the problem of enhanced gravity in spiral 
galaxies with the Vaidya metric changes to the problem of the energy source 
and the physical nature of its emission. Physically, it is very different from 
the Newtonian solution, which is the nonluminous matter. 

Presently it is a standard assumption that galactic dynamics is governed 
by dark matter. To prove the dark matter hypothesis a number of experi- 
ments start to search for neutralino, the best supersymmetry candidate for 
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WIMP. Also, astrophysical observations of dark matter in elliptical galaxies 
have been attempted by PN.S collaboration [2] with planetary nebulae as 
a tracer of gravity Surprisingly gravity of those galaxies is adequately de- 
scribed by luminous matter only , a result described as a "missing missing 
mass" problem [2]. If the results obtained by PN.S collaboration are correct 
than dark matter in elliptical galaxies is at least differently distributed than 
in spiral galaxies, with only trace amount inside inner 5-6 effective radii. 

A radical proposal is the Milgrom's modified Newtonian dynamics (MOND) 
hypothesis, which postulates new gravitation law for very weak acceler- 
ations. This proposal, employed as a phenomenological model, is quite 
successful in explaining spiral galaxies dynamics. MOND also explains the 
dearth of dark matter in elliptical galaxies [8] observed by Romanowsky et 
al. [2]. The radial dependence of the MOND acceleration is the same as in 
non-Newtonian acceleration (4.14) for the Vaidya metric. 

One can notice that the Vaidya metric is an example of metric considered 
recently by Lake [9] that can give flat rotation curves of spiral galaxies. 

We have shown here that magnetic fields in spiral galaxies can play 
a crucial role in determining the rotation curves. The solution (3.6) shows 
that magnetic field which fully accounts for a flat rotation curve has toroidal 
component compatible with observed magnetic fields in spiral galaxies. It 
may not be a good approximation to completely suppress magnetic field in- 
fluence when studying the physical origin of flat rotation curves. In realistic 
description observed magnetic field influence should be subracted before fit- 
ting gravitational potential generated by assumed dark matter halo. The 
importance of magnetic field contribution to flat rotational curves of spiral 
galaxies has been recently discussed in Ref.[10]. 

It is also worth to notice that pure magnetic mechanism of flat rotation 
curves in spiral galaxies could explain simultaneousely why these curves in 
elliptical galaxies are Keplerian. It is because there are no regular magnetic 
fields in elliptical galaxies. 
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